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Abstract 

In this work we perform gfobal fits of microscopic decoherence models of neutrinos to all avail- 



' able current data, including LSND and KamLAND spectral distortion results. In previous works 

00 ' 

CN ■ on related issues the models used were supposed to explain LSND results by means of quantum 

O , 

. gravity induced decoherence. However those models were purely phenomenological without any 



underlying microscopic basis. It is one of the main purposes of this article to use detailed mi- 
croscopic decoherence models with complete positivity, to fit the data. The decoherence in these 
• models has contributions not only from stochastic quantum gravity vacua operating as a medium, 

r~| , but also from conventional uncertainties in the energy of the (anti) neutrino beam. All these con- 

^ tributions lead to oscillation-length independent damping factors modulating the oscillatory terms 

}^ • from which one obtains an excellent fit to all available neutrino data, including LSND and Kam- 

land spectral distortion. The fit is much superior to all earlier ones. It appears that the results of 
the fit are most naturally interpreted as corresponding to conventional energy uncertainties. This 
represents a radical departure from previous analyses where the neutrino data (including LSND 
but not Kamland spectral distortion) were regarded as evidence for quantum gravity decoherence. 

PACS numbers: 14.60.Pq, 04.60.-m 
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I. INTRODUCTION 



The theory of Quantum Gravity (QG) is still elusive. In some theoretical models, the 
phenomenon of space-time 'foam', invoked by J. A. Wheeler may be in place; according 
to this picture the singular microscopic fluctuations of the metric, give the ground state of 
QG the structure of a 'stochastic medium'. The medium has the profound effect of leading 
to decoherence of quantum matter as it propagates. This may have experimentally observable 
consequences in principle 

One of the basic effects of decoherence is the presence of damping factors in front of 
the oscillatory terms. However, one should be very careful when interpreting decoherence 
effects, if observed in an experiment, because ordinary matter can easily 'fake' decoherence 
effects, especially the damping exponents 0,0]. For instance, uncertainties in the energy of 
a neutrino beam P| , which are associated with ordinary physics, and have nothing to do with 
'fuzziness' of space time, do reproduce a damping exponent similar to that encountered in 
Lindblad decoherence models p|. Of course, stochastic quantum gravity effects can induce 
such uncertainties in the energy beam, and hence contribute to the damping exponent 
themselves but such effects are usually subleading. Thus, one should know the energy 
of the beam with high precision in order to eliminate 'fake' decoherence effects and probe 
quantum gravity effects sufficiently well. 

In ref. H, henceforth referred to as I, we have attempted to fit the available neutrino 



data, including LSND results [2I|, using phenomenological decoherent models with mixing in 
all three generations of neutrinos. Such fits extended earlier similar attempts to study deco- 
herence with two-generation neutrino models S|. In I it was seriously entertained that the 
decoherence might be attributable to environmental entanglement with the quantum gravity 
foamy vacuum, and could be distinguished from ordinary-matter-induced decoherence P|. 

A simplified model, of Lindblad type has been used for the fit, following earlier work 
in P]. The model of |9] involved a phenomenological diagonal decoherence matrix, and in 
the fit of I, decoherence was assumed to be dominant only in the antineutrino sector. This 
assumption was made in order to fit the LSND results [7| pointing to significant ^ Vf. 
oscillations, but no significant evidence of oscillations in the particle sector. In this way 
a fit was made to a three generation model with the LSND "anomalous" result, without 
introducing a sterile neutrino. The possibility of strong CPT violation in the decoherence 
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sector, allowed for an equality of neutrino and antineutrino mass differences in agreement 
with atmospheric and solar neutrino data. 

The particular choice of j^, which yielded an extremely good fit to all available neutrino 
data involved mixed energy dependence for the (antineutrino-sector) decoherence coeffi- 
cients, some of which were proportional to the neutrino energies E, while the rest were 
inversely proportional to it, oc 1/E. In I, the coefficients proportional to 1/E were inter- 
preted as describing ordinary matter effects, whilst those proportional to E were assumed 
to correspond to genuine quantum-gravity effects. The latter increase with the energy of 
the the (anti)neutrino is consistent with a larger back reaction effect on quantum space-time 
and hence with a larger decoherence. 

The strong difference assumed in I between the decoherence coefficients of the particle and 
antiparticle sectors, although not incompatible with a breakdown of CPT at a fundamental 
level P], appears at first sight somewhat curious, and in fact is unlike any other case of de- 
coherence in other sensitive particle probes, like neutral mesons, examined in the past jlfll |. 
There, the oscillations between particle and antiparticle sectors, necessitate a common de- 
coherence environment between mesons and antimesons. If one accepts the Universality 
of gravity, then, the sample point of I seems incompatible with this property. Moreover, 
there are two more problematic points of the fit in I, which were already discussed in that 
reference. The first point concerns the complete positivity of the model. In the ad hoc 
diagonal form of the decoherence matrix, used in I, was postulated , without a discussion of 
the necessary conditions required in the Lindblad approach to guarantee complete positiv- 
ity. Indeed the particular choice of the decoherence parameters of I, did not lead to positive 
definite probabilities for the entire regime of the parameter space of the model, although 
the probabilities were positive definite for the portion of the parameter space appropriate 
for the various neutrino experiments used for the fit. Specifically, it was found in I, that 
with the particular choice of the decoherence parameters in the (antineutrino sector), one 
obtains positive-definite transition probabilities for energies restricted to E > 0{1 MeV). 
The second, and more important point, is that the choices of decoherence parameters of I 
were good for all the neutrino experiments available at the time, but unfortunately it could 
not reproduce the spectral distortion observed by the KamLAND experiment JAJ, whose 
first results came out simultaneously with the results of I. 

The aim of this article is two fold. One is to rectify the above points, and present a 
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novel fit, using as in I the simplified three-generation Lindblad model of decoherence of 
0], but crucially amended so as to respect the general conditions among the coefficients 
necessary to guarantee complete positivity in the entire parameter space. We shall show 
below, that it is possible to find such a consistent Lindblad model consistently, for which 
the fit to all available neutrino data is excellent, including the spectral distortion seen by 
KamLAND, and the LSND results j^. This substantially extended decoherence model (in 
comparison to I), constitutes therefore the first mathematically consistent three-generation 
neutrino decoherence model of Lindblad type which fits all the available data, including 
spectral distortion seen by KamLAND. 

The second and more significant aim was to give a microscopic and physically motivated 
model which would fit into the general scheme of the linear Lindblad decoherence. In 
this way the constraints obtained from the phenomenological fit can be examined to check 
consistency with values that can be deemed reasonable for phenomena originating from 
quantum gravity. Such a comparison has not been done before and leads to a major shift in 
our views concerning decoherence due to quantum gravity. In fact, as we shall argue below, 
the most natural explanation of the fit seems to be provided by energy uncertainties in the 
(anti)neutrino beam, due to conventional physics. Several microscopic quantum space-time 
models, that we have examined in this work, yield too small effects to reproduce the result 
of the fit. 

In the present work, the decoherence parameters in the model are assumed to be the same 
in both neutrino and antineutrino sectors, consistent with the above-mentioned universal 
property of a quantum-gravity environment. In this sense, we assume that the LSND result 
is correct in both channels, although their observed excess of z/g events is not corroborated (at 
the same level at least) in the neutrino channel. This is to be contrasted with the approach 



of I, where following 



12j |. only the evidence in the antineutrino sector was considered. 



The structure of this article is the following: in section 2 we review the basic theory 
of Lindblad decoherence, and specify the conditions for complete positivity in the type of 
model for decoherence used in 0| and in I. We examine the limitations imposed on the 
parameter space of the model in I in order to guarantee complete positivity, and then we 
construct a modified model, in which one obtains positive definite transition probabilities for 
the entire regime of the parameter space. The decoherence implies exponential damping with 
time (oscillation length), which violates microscopic time irreversibility, irrespective of CP 



4 



ation. In order to agree with the experimental 
111 we require such exponential damping factors 



properties of neutrinos, and hence CPT vio 
results of KamLAND on spectral distortions 
to imply a modulation in front of some of the oscillatory terms giving rise to a modification 
in the (survival) transition probabilities of order per mil. We obtain in section 3 stringent 
constraints on the exponents of the damping factors. The sample point that fits all available 
data, including LSND,is discussed in detail in section 3. This model-point corresponds to an 
exponent of the decoherent damping factors which is independent of the oscillation length. 
An attempt to explain such a result in terms of microscopic models of stochastic space time 
foam is given in section 4. However, as we show there, explanations based on conventional 
physics such as the uncertainties in the (anti) neutrino beam energy, are definitely much 
more plausible. The present data when interpreted in terms of a microscopic model make 
quantum gravity an unlikely candidate for the origin of decoherence (claimed to be observed 
in our fits). Finally, conclusions and outlook are presented in section 5. 

II. LINDBLAD DECOHERENCE AND TRANSITION PROBABILITIES: A RE- 
VIEW 

In this section we present the details of the calculation for transition probabilities of 
three generations of neutrinos, where complete positivity is maintained within the (linear) 
Lindblad approach In this framework the general evolution equation of the p density 
matrix, representing a (spinless) neutrino state reads : 

dtp = L[p] (2.1) 

where there are conditions on the decoherence contribution to L which guarantee complete 
positivity of the probabilities as they evolve in time. The spin of the neutrino will not play 
an important role in constraining the decoherence sector by comparing with experimental 
data, and hence we shall present a formalism based on scalar particles. Detailed studies of 
Dirac and spinless neutrinos have been performed in ; as explained there the inclusion 
of spin does not affect qualitatively the main decoherence effects which is the damping of 
oscillation probabilities. In section 4, where we attempt to interpret the time (i.e. oscillation 
length) dependence and order of magnitude of the decoherence parameters, we shall present 
a more detailed discussion on the results of 
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With the above in mind, we commence our analysis with a theorem due to Gorini, Kos- 
sakaowski and Sudarshan lljl on the structure of L, the generator of a quantum dynamical 

nn 

semi-group |5|,|1J]. For a non- negative matrix Cki (i.e. a matrix with non- negative eigenval- 
ues) such a generator is given by 



k,l 



where H = if^ is a hermitian Hamiltonian, {Fk,k = 0, 



(2.2) 

,77,^ — 1} is a basis in M„(C) 



such that Ft 







Tr(Ffc) = OWk ^ and Tr(i^"^F_,) = 6ij. In our application we can 



take Fj 



2 



(where Aj are the Gell-Mann matrices) and satisfy the Lie algebra [Fi,Fj 



i^^j. fijkFk, (yi = 1,...8) fijk being the standard structure constants, antisymmetric in all 
indices. 

Without a microscopic model, in the three generation case, the precise physical signifi- 
cance of the decoherence matrix cannot be fully understood. In this work we shall consider 
the simplified case in which the matrix C = {cki) is assumed to be of the form 
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(2.3) 



As stated above, positivity can be guaranteed if and only if the matrix C is positive and 



hence has non-negative eigenva! 



simplification has been used in 



ues. We have also taken C to be symmetric. A similar 



and in I, to yield an economic decoherence model which 



can be compared with experimental data. However, as discussed in detail in such special 
choices can be realised for models of the propagation of neutrinos in models of stochastically 
fluctuating environments where the decoherence term corresponds to an appropriate 
double commutator involving operators that entangle with the environment. The quantum- 
gravity space time foam may in principle behave as one such stochastic environment 

m 



13l |. and it is this point of view, that we will critically examine in this work. In section 4, we 
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shall discuss the viability of the interpretation of the fit in terms of such microscopic models 
of space-time foam. 

Since the are Hermitian, we can rewrite the expression for L[p\ as 

L[p\ = -z[H, p] + V[p] = -z[H, p] + \Y1 ^'^^ (t^'^'"' + t^-^' Z'^^]) (2.4) 

k,l 

After standard manipulations, we may write the non-Hamiltonian decoherence part as 
V[p] = i ([F,, [p,Fi]] + {F,, [p,Fi]} - [Fi, [F,,p]] + {Fi[F,, p]}) + ^{p, [F,, F,]} (2.5) 
On using the expansion p = "^.piFi, this expression can be written 

+ifkim{-5lm. + ^flmjFj) + i2fklm{-Sim + ^^imjFj)^ (2.6) 

We note that the only terms which contribute are {—fumfkmj + fkimfimj) Fj. 

We follow the basic notation and express the time evolution of the density matrix 

as 

Pk = ^ifijk + '^kj)pj = ^ Mkjpj (2.7) 

j i j 

where we have 

'^ij ^ ^ ^ ( filmfkmj ~l~ fkimfimj^ ■ l^'^) 
k,l,Tn 

Using the values of the structure constants fijk of the SU(3) group, appropriate to the three 
generation case being examined here, we arrive at: 

_ 1 / 1 

T^U — ~2 [ '^22 + '^33 + ^('^44 + C55 + Cge + 077^ 

T^22 = ~2 I 4 ('^'^^ "^^^ 

_ 1 / 1 

^33 — ^2 ( '^11 '"22 + ^^'^44 + C55 + C66 + C77, 

_ 1 / 1 Q ^ 

^44 — ~2 1 '"^^ ~'~ 4''""'^^ ~'~ '"^^ ~'~ ""^^ ~'~ ""^^ ~'~ '"^'^ ~'~ '"^'^^ — 2~'"'^^ 

^^55 = "2 I C44 + ^(Cll + C22 + C33 + C66 + C77 + 3C88) + -^^33 



2^83 



III, Q A 

C77 + ^I'^ll + '^22 + C33 + C44 + C55 + oC88j 2~'''^^ 

1(1, o \ ^ 



--(C44 + C55 + C66 + C77) 



V 



38 



^(C44 + C55 - C66 - C77) 



(2.9) 



or conversely, 

Cll 

C22 
C33 
C44 

C55 

C77 



3^^88 + 
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11 



^1)88 
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11 



1^22 "~ ^^33 
+ - ^^33 

1^22 + ^^33 



^^55 + 



44 



3 ^ 



1^55 - 1^44 - ^V3^^38 - ^^^8S 
-P77 + 'T>m + ^^^38 - ^1^8 

2 2 
D77 — Pee + ^2 "^38 ~ 3^88 



2^ 2^ 2^ 

— 1^55 - -^^77 - -I'e 
3 3 3 



^2^44 + 



^Pn + \V22 + ^1^33 



C38 



(2.10) 



The simplified form of the Cj^ matrix given in ()2.3p implies a matrix Ljj of the form: 



^1^11 
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(2.11) 



where we have used the notation A,' 



2p 



The corresponding eigenvalues are: 
1 



Ai 
A2 
A3 
A4 
A5 
As 
A7 
A« 
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2 



[(^^44 + v. 
[(^^44 + 
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2 

^[(©66+^^77 
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88 
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22 



33 



+ m 



38J 



44 - I?55)2 - 4A? 
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4A? 



13J 



^77 J 



4Ai3] 



(Pee - Vrry 



4A^3] 



(P33 - ^^88) 



4^^i8] 



(Vn + 1?22 
(Pn + ^^22 
(Dn + ^?22 

(P44 + ^^55 
(1^44 + 1^55 
(©66 + P77 
(©66 + ^^77 
{Vu + V22 



-^12 
+ ^12 
— ^38 
-^^13 
+ ^13 
-^23 
+ ^23 
+ ^38 



(2.12) 



The probabihty of a neutrino of flavor z/q,, created at time t = 0, being converted to a flavor 



Up at a later time t, is calculated in the Lindblad framework 

1 1 



to be 



P..^.,(t) = Tr[p-(t)p/5] = - + -$^e^'^-*A./^^Vj(0)pf . 



where the matrix D and its inverse are 



/ Aj— X'22 A2— X'22 
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(2.13) 



(2.14) 
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and 
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(2.15) 



/ 

It will be sufficient to look explicitly at the k = 1 and k = 2 terms in the sum of the right 
hand side of ()2.13|) . since by block symmetry the other terms will be of the same form. 
For the k = 1 term we have: 



e 2 e 2 



a /3 + ^22 + ^12 , a /3 ~^22 + ^^11 + ^12 
PlPl +P2P2 
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, a /3A12 I a 13 

"12 "12 



(2.16) 



Likewise for A; = 2 we have 

e"^* (A2/^2VPi (0)pf ) = e^^* [p-lpfDuD^,' + p^pf ^22/^22' + P2PfDi2D^i + P?pf I^22/^2"i') 



(P11+P22)* "12* 
e 2 e 2 



■ „ f} -V22 + Vii + ni2 ^ „ /3 (-1^22 + 1^11 - ^^12) 
PlPl ^7^ + P2P2 
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12 



, a a~^12 , a B^12 
+ P2P1-P^ + P1P2 



n 



12 
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Upon combining equations ()2.17|) and ()2.18|) we obtain: 



a /3 I a 13 
P1P1+P2P2 



Q-^r2^ _|_ gf^i2| 



2Ai2(p?pf - P^pD + AP21 (P?P? - Pfp^) \ / e-^'-5 - e 



,/5 .a 



(2.17) 



(2.18) 



+ 



12 



10 



As mentioned earlier, by block symmetry we can see that the other terms will be of the 
same form. 

We thus obtain for the relevant probability: 



1 1 
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P4Pf + P5 P5 



12 
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13 
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'2A23(p?P? - P?pS + AVre [p^pI - p^f^) 
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^23 



2D38(P?Ph - /^?/'?) + ^2^83 ( P?P3 - f>"Pi 



,0 



38 



g-ni2| _ £^12 1 



^(r'ii+i?22)| 



,(r'44+I'55)| 



g-n23| _ Q^23^ 



3 --.,N o _ ,,--.,8 2 



,(I?33+f88)| 



Above we have used the notation that AP^j = Da — Vjj. We have assumed that 
2|Ajj| < lADjjl with the consequence that Q^j, ij — 12,13,23 is imaginary. However, 



Q38 = v^(X'33 - VssY + 4X'|8 will be real. Thus, the final expression for the probability 
reads 



1 1 
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3 2 



„ p I l^i2|A ^ / APsipfp? 
PiPfcos(^J + (^^^^^ 
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11 



^ ' Vl^s ^^^^ 



,(©33+1)88)1 



(2.19) 



On using the relations 



Po 



2 

p° = 2Re{U:,U^2) 



p° = -2/m(f/:if/,2) 

Ps = — |f^Q2P 

p° = 2Re{U:,Ua3) 
pt = -2Im{U:,Ua3) 
p° = 2ReiU:,U^s) 



= -2Im{U:,U, 



a3) 

Ps = )jl{\Ual\^+\Ua2\^-2\U^s\^). 

we can readily see that for a real U matrix (i.e. no CP violating phases) the relevant 
probabilities are bounded. Indeed, there is no danger of the cosh and sinh terms blowing 
up with time t, as we always have Q^s < P33 + X^gg. We can see this by checking that 

(D33-^^88)'+4Pi8 < Pi3+Pig + 2^332^88 , -2^33^^88 + 42^18 < 22^33^^88 , and Pfg < 1^332^88, 

which are automatically satisfied, as can be readily checked from the relevant expressions 
(EH). 

We are now ready to discuss the fit to the experimental data. This is done in the next 
section. 



III. FITTING THE EXPERIMENTAL DATA 

In order to check the viability of our simplified scenario, we have performed a com- 
parison (as opposed to a fif ) fo SuperKamiokande sub-GeV and multi GeV data (the 40 
data points that are shown in Figure 1), CHOOZ data (15 data points), KamLAND (13 data 
points, shown in Figure 2) and LSND (1 datum), for a sample point in the vast parameter 
space of our extremely simplified version of decoherence models. Rather than performing 
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a x^"fit (understood as a run over all the parameter space to find the global minimum of 
the function) we have selected (by "eye" and not by x) ^ point which is not optimised 
to give the best fit to the existing data. Instead, our sample point must be regarded as a 
local minimum around a starting point chosen by an educated guess. It follows then that 
it may be quite possible to find a better fitting point through a complete (and highly time 
consuming) scan over the whole parameter space. 

To simplify the analysis and gain intuition concerning the rather cumbersome expressions 
for the transition probabilities, we have imposed. 





= V22 , 


P44 


= 2^55 






Vi3 


= 1^88 




^^38 = ^^83 


= . 





(3.1) 

implying a diagonal "D-matrix. 

Later on we shall set some of the Va to zero. Furthermore, we have also set the CP 
violating phase of the KMS matrix to zero, so that all the mixing matrix elements become 
real. 

With these assumptions, the complicated expression for the transition probability (|2.19p 
simplifies to: 



+ p2p'. cos e--* + p-pI cos e^-* 

+ (p3P^ + P8Pf)cosh(%^)e-33*| (3.2) 

for both, neutrino and antineutrino sectors. 

As indicated by the state-of-the-art analysis, masses and mixing angles are selected to 
have the values 

Amj^ = 7 ■ 10-5 eV2, 
Am^g = 2.5 ■ 10-3 eV2, 
^23 = vr/4, ^12 = .45, = .05 



For the decoherence parameters we find 

T>33 = 2^66 = , Vu = V22 = Vu = 1^55 = - 

13 



1.3-10-2 



(3.3) 



in units of 1/km with L = t the oscillation length. The 1/L-behaviour of Pn, implies 
oscillation-length independent Lindblad exponents. We shall attempt to interpret this result 
in the next section. 

The complete positivity of the case defined by (j3.1|) and (j3.3|) is guaranteed; this follows 
from the fact that the solutions for Cij in terms of T>ik (cf. (|2.1H) ) is such that the only 
non-zero elements are cgs = cas/v^- Such a C-matrix has only non- negative eigenvalues 

C - matrix eigenvalues = (0, 0, 0, 0, 0, 0, 0, -8 Pii/3) (3.4) 

In summary we have introduced only one new parameter, a new degree of freedom, by 
means of which we shall try to explain all the available experimental data. It is important to 
stress that the inclusion of one new degree of freedom by itself does not guarantee that all the 
experimental observations can be accounted for. Indeed for situations without decoherence 
the addition of a sterile neutrino (which introduces four new degrees of freedom -excluding 
CP violating phases) seemed to be insufficient for matching all available experimental data, 
at least in CPT conserving situations jioj l. 

In order to test our model with this decoherence parameter, we have calculated the zenith 
angle dependence of the ratio "observed-events/ (expected-events in the no oscillation case)" , 
for muon and electron atmospheric neutrinos, for the sub-GeV and multi-GeV energy ranges, 
when mixing is taken into account. Since matter effects are important for atmospheric 
neutrinos, we have implemented them through a two-shell model, where the density in the 
mantle (core) is taken to be roughly 3.35 (8.44) gr/cm^, and the core radius is taken to be 
2887 km. 

The results are shown in Fig. 1, where, for the sake of comparison, we have also included 
the experimental data. As can be easily seen the agreement is remarkable. 

As bare eye comparisons can be misleading, we have also calculated the value for each 
of the cases, defining the atmospheric as 

10 / pexp _ pth \2 
M,S a=e,ii i=l 

Here the statistical errors, the ratios R^^i between the observed and predicted signal 

can be written as 

KJ = K7/<f (3-6) 
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FIG. 1: Decoherence fit. The dots correspond to SK data. 

(with a indicating the lepton flavor and i counting the different bins, ten in total) and M, S 
stand for the multi-GeV and sub-GeV data respectively. For the CHOOZ experiment we 
used the 15 data points with their statistical errors, where in each bin wc averaged the 
probability over energy. For the KamLAND experiment, their 13 data points have been 
used for a flxed distance of Lq = 180 km, as if all antineutrinos detected in KamLAND were 
due to a single reactor at this distance and plotted in figure 2 while for LSND one datum 
has been included. 

The results are summarised in Table 1, where we present the comparison for the model 
in question and the standard scenario (calculated with the same program). 

From the table it becomes clear that our simplified version of decoherence in both neutrino 
and antineutrino sectors can easily account for all the available experimental information, 
including LSND. It is important to stress once more that our sample point was not obtained 
through a scan over all the parameter space, but by an educated guess, and therefore plenty 
of room is left for improvements. On the other hand, for the mixing-only/no-decoherence 
scenario, we have taken the best fit values of the state of the art analysis and therefore no 
significant improvements are expected. 

As we have seen, the decoherence effects suffered by our model, are just an overall sup- 
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FIG. 2: Ratio of the observed spectrum to the expectation versus Lq/E for our decoherence 
model. The dots correspond to KamLAND data 





decoherence 


standard scenario 


SK sub-GeV 


38.0 


38.2 


SK Multi-GeV 


11.7 


11.2 


Chooz 


4.5 


4.5 


KamLAND 


16.7 


16.6 


LSND 


0. 


6.8 


TOTAL 


70.9 


77.3 



TABLE I: obtained for our model and the one obtained in the standard scenario for the different 
experiments calculated with the same program. 

pression on some of the oscillatory terms modifying the transition/survival probabilities at 
the per mil level (to account precisely for LSND, a per mil evidence) and therefore, no effect 
is expected (or found) in the oscillation dominated physics, where the level of precision is at 
the percent level, at most. We are guaranteed then to have an excellent agreement with solar 
data, as long as we keep the relevant mass difference and mixing angle within the LMA-I 
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region, something which we certainly did. Thus, there is no need to include these data on 
our fit. 

At this point, a word of warning is in order. Although from the table, it seems that the 
decoherence model we are presenting here and the standard no-decoherence scenario provide 
equally good a fit, i.e. while the former has a x^/DOF = 70.9/63 the latter has a x^/DOF 
= 77.3/64, both quite "acceptable" from the statistical point of view, one must remember 
that only the decoherence model can explain the LSND result. This fact, however, gets 
blurred in the total because LSND is represented by only one experimental point with a 
poor precision. 

Before closing this section, it is worth revisiting the models of I, in order to understand 
in the above context, the failure of complete positivity in certain regions of the parameter 
space. In that case, the following restrictions on the decoherence matrix (which was also 
diagonal, as in the case p.ip above) had been imposed j^: 

■E = -A, 

-24 



1^11 


= V22 


= Vu 


= 2^55 = 


-2 ■ 




= V-j-j 


= 2^33 


= ^^88 = 


-10 




= Vs3 


= 0, 







10-^7^ = -B , 

(3.7) 

leading to a solution for the c-matrix c^s = | V^A — | y/SB, C55 = | i?, C44 = | i?, Cgg = 



^ A, cqq = 1-6, C22 = iB, C33 = -A/SB + 2 A, cjj = f 5, Cn = | 5 (all other matrix 

8/1 2d2d2d2d2d2d2 



3 ^^5 '-'DD 3 ' 3 ' ■J'J j-/^-^ I ^j-i, 3 ' 3 

entries zero) such that the pertinent eigenvalues (—2 -B+| ^4, | i?, | i?, | i?, | i?, | i?, | 5, | B) 
which are not positive for arbitrary values of A and B. The positivity condition can be 
obtained by demanding positivity of the first eigenvalue, i.e. 2 ■ lO"^^/^' < y ■ 10~^^ ■ E, 
where E is in units of GeV, which leads to the condition E > 0{1 MeV) which was the 
condition found in I. 



IV. ATTEMPT AT INTERPRETING THE FIT 



The microscopic origin of the "observed" decoherence effects, according to our fit above 
( ()3.3|l . ()3.4j) ). may not be unique. In fact there can be many contributions to the decoherence- 
induced (oscillation-independent) damping ()3.4j) . which modulates the oscillatory terms, 
arising from a variety of effects, ranging from microscopic quantum-space time fluctuations 
('stochastic quantum gravity foam' jl,]), to ordinary matter effects, e.g. uncertainties in the 
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energy and oscillation length of the (anti)neutrino beam. It is the purpose of this section 
to attempt and disentangle these very different in nature contributions, and in particular to 
estimate their plausible order in terms of microscopic theoretical models and see which one 
gives the order specified by the fit. 

To understand the results (|3.3|) . (j3.4|) . in connection with either stochastically fluctuat- 
ing quantum-gravity space-time-foam models , or energy-uncertainty driven deco- 
herence jsj, it suffices to restrict our discussion to the simplified (but phenomenologi cally 
realistic ) case of three neutrino families, but with dominant mixing only between 12, 23 |l3[, 
with mixing angles 612 = 623 = and ^^13 = 0. The three generation case with full mixing 
does not affect qualitatively the form of the damping exponents used to fit of the oscillation 
experiments, and hence we are free to use results on theoretical models from this simpler 
case, in order to interpret qualitatively the above result. 

A. Stochastic Quantum-Gravity Models 

We commence our analysis with models of neutrinos propagating in a quantum-gravity 
ground state. To this end, consider the propagation of such a neutrino system, in a medium 
whose density stochastically fluctuates. For our purposes the medium is taken to be a quan- 
tum space time foam with fluctuating densities of charged black-hole/anti-black-hole 
pairs produced by the vacuum and being absorbed by it within Planckian time scales 12 1. 
Such a case, will not produce any vacuum charge on average, but the associated density 
fluctuations will produce vacuum fluctuations in electron currents with which electron neu- 
trinos will interact. An inherent CPT violation may result in asymmetries between particle 
antiparticle sectors, as far as the appropriate interactions of the (anti) neutrinos with these 
currents are concerned. 

Such asymmetries can produce a bias in flavour of, say, the electron current fluctuations 
in the space-time foam vacuum, with the result that a "gravitationally-induced" MSW 
effect is in place , with small contributions to the standard oscillations due to bare mass 
differences between neutrinos, that could be due to non-gravitational physics 0. The 
precise reason why a preference to electron current fluctuations as opposed to positron ones 
cannot be answered at this stage, since a microscopic model of quantum space time foam 
is lacking. However, the situation is not incompatible with the intrinsic CPT violation (to 



18 



be precise microscopic time irreversibility, unrelated to CP properties, according to which 
the generator of time reversal operations is ill-defined) characterising such problems, where 
a proper scattering matrix cannot be defined jl^ . 

In such a case, the evolution equation of the density matrix p of the neutrino probe in- 
volves a time-reversal (CPT) breaking decoherence matrix of a double commutator form jisj ], 

dtp = t[p,H,ff]-n'[Hj,[Hi,p]] (4.1) 

where {n{r)n{r')) = Q'^nl6{r — r') denote the stochastic (Gaussian) fiuctuations of the 
density of the medium and Hp f f = Ho + Hj, Hq being the standard Hamiltonian, and Hi 
an MS W-type interaction 3, Q] • This double-commutator decoherence is a specific case 
of Lindblad evolution, of the type considered in previous sections, with a C-matrix of the 
form: 



C 



(4.2) 



[ hj hihs hihs \ 

0000 

hih3 hj hhs 

0000 

0000 

0000 

0000 

\hihs hshs hj J 

with hi = {tty^ - a^^)sm{2e), = {a^^ - a^,^) cos{26), and h, 
indeed has positive eigenvalues for real hi. The difference - 
average density of the medium no- 

We note at this stage that, for gravitationally-induced MSW effects (due to, say, black- 
hole foam models as in ^iZ]), one may write 

Atte^ = a^, - a^,^ oc Gntiq (4.3) 

with Gn = 1/Mp, Mp ~ 10^^ GeV, the four-dimensional Planck scale. This gravitational 
coupling replaces the weak interaction Fermi coupling constant Gp in the conventional MSW 
effect. This is the case we shall be interested in for the purposes of this work. 

In such a case the density fluctuations are therefore assumed small compared to other 
quantities present in the formulae, and an expansion to leading order in fl"^ is in place. 



^^^^^^ . This matrix 

v3 

is proportional to the 
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Following then the standard analysis, outlined above, one obtains the following expression 
for the neutrino transition probability z/g ^ z/^ in this case, to leading order in the small 
parameter <^ 1: 



12 I 4p5/2 p3/2 



a2 A2 



_g-Aa2^n2,(i+_i^(,„3(4,)„l)) ^Qg2(^) Sin2(2^)- 



12 



2r 

-e -r oos\ey-^ ' + 2 cos'W (4-4) 

where T = (Aoe^ cos(2^) + Au? + Aal^ sin\2d) , Au = ^ ■ 

From ()4.4j) we easily conclude that the exponents of the damping factors due to the 
stochastic-medium-induced decoherence, are therefore of the generic form, for t = L, the 
oscillation length (in units of c = 1): 

A 2 ..2 A (C0S(4^) - 1)\ 

exponent ~ Aa^nh i^l + ' '-j (4.5) 

The reader should note at this stage that, in the limit A12 0, which could characterise the 
situation in j^, where the space-time foam effects on the induced neutrino mass difference are 
the dominant ones, the damping factor is of the form exponentgj.a^yjtj^tionai msw ~ 
with the precise value of the mixing angle 9 not affecting the leading order of the various 
exponents. However, in that C3iSG, clS follows from ()4.4|) . the overall oscillation probability 
is suppressed by factors proportional to A^g, and, hence, the stochastic gravitational MSW 
effect 5 although in principle capable of inducing mass differences for neutrinos, however 
does not suffice to produce the bulk of the oscillation probability, which is thus attributed 
to conventional flavour physics. 

In what follows, therefore, we assume the case where A12 ^ Aa^.^, and this is the case we 
shall compare with the results of our experimental fit above. The result of the fit ()3.3|) . ()3.4p . 
then, implies that the above decoherence-induced damping exponent ()4.5|) is independent of 
L and actually we have, to leading order in Aae^/Ai2 ^ 1 (re-instating dimensions of h, c): 

n^(Aae^)2 (^1 + ^Q<4g)-1 ^ _ ^ ^ 2.56 X IQ-^^ GeV ■ km. (4.6) 

This in turn implies that in this specific model of foam, the density fluctuations of the 
space-time charged black holes is such that, for maximal mixing, say, 9 = n/A assumed 
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for concreteness, and for L ~ 180 Km, as appropriate for the KamLAND experiment, the 
decoherence damping factor is P = Vt^G\n^ ~ 2.84 ■ 10^^^ GeV, if the result of the fit 
is due exclusively to this effect (note that the mixing angle part does not affect the order 
of the exponent). Smaller values are found for longer L, such as in atmospheric neutrino 
experiments. 

The independence of the relevant damping exponent from the oscillation length, then, 
implied by our fit above, may be understood as follows in this context: In the spirit of 17|, 
the quantity GnHq = where ^ ^ 1 parametrises the contributions of the foam to the 

induced neutrino mass differences, according to our discussion above. Hence, the damping 
exponent becomes in this case ^^fi^(Am^)^ • L/E"^. Thus, for oscillation lengths L we have 
L^^ ~ Am'^/E, and one is left with the following estimate for the dimensionless quantity 
^'^Am'^Q'^/E ~ 1.3 ■ 10~^. This implies that the quantity Q"^ is proportional to the probe 
energy E. In principle, this is not an unreasonable result, and it is in the spirit of jl7|. 
since back reaction effects onto space time, which affect the stochastic fluctuations Q"^, are 
expected to increase with the probe energy E. However, due to the smallness of the quantity 
Am'^/E, for energies of the order of GeV, and Am^ ~ 10^'^ eV^, we conclude (taking into 
account that ,^ <^ 1) that Q'^ in this case would be unrealistically large for a quantum-gravity 
effect in the model. 

We remark at this point that, in such a model, we can in principle bound independently 
the Q and Uq parameters by looking at the modifications induced by the medium in the 
arguments of the oscillatory functions of the probability ()4.4|1 . that is the period of oscillation. 
Unfortunately this is too small to be detected in the above example, for which hi <^ A12. 

The result of the fit, however, may be interpreted more generally, as implying independent 
of the oscillation length L exponents in the decoherence exponential suppression factors 
in front of the oscillatory terms in the transition probabilities. In this sense, the bound 
fl3.3|) . ()3.4|) determined by the fit above, can be applied to other stochastic decoherence 
models, for instance the one discussed in in which one averages over random (Gaussian) 
fluctuations of the background space-time metric over which the neutrino propagates. 

In such an approach, one considers merely the Hamiltonian of the neutrino in a stochas- 
tic metric background. This is one contribution to decoherence, since other possible non- 
Hamiltonian (like the Lindblad terms above) interactions of the neutrino with the foam are 
ignored. In this case, one obtains transition probabilities with exponential damping factors 
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in front of the oscillatory terms, but now the scaling with the oscillation length (time) is 
quadratic consistent with time reversal invariance of the neutrino Hamiltonian. For in- 
stance, for the two generation case, which suffices for our qualitative purposes in this work, 
we have: 

^g»K-c.2)t^ = k e H H ))x 

!!_^(^^,+^,+^3+^,)+ ^'--^'^(-^-i) (12^,+2a2-2a3)V^ 



(4.7) 

where k is the neutrino energy, cXj ,i = 1,...4 parameterise appropriately the stochastic 
fluctuations of the metric in the model of Il3|, T = — r-^, |T| <^ 1, and k"^ ^ m? , mn, 
and 

2+ = i (m^ + T(l + cos2e){ml - m^) + T^{m\ - ml) sin^ 2^) 



1 







- + T(l - cos2^)(m? -m^) - T2(m? -m2)sin2 2^) . (4.8) 



le oscillation period, as well 
1 

13| . The latter is attributed 



2 

Note that the metric fiuctuations-cij induced modifications of t^ 

as exponential e"'-'' -'*^ time-reversal invariant damping factors 

to the fact that in this approach, only the Hamiltonian terms are taken into account (in a 

stochastically fiuctuating metric background), and as such time reversal invariance t ^ —t 

is not broken explicitly. But there is of course decoherence, and the associated damping. 

We, then, observe that the result of the fit above, ()H.Hj) . ()H.4|l . implying L-independent 

exponents in the associated damping factors due to decoherence, may also apply to this 

case, implying for the damping exponent: 

'm?— mn)^, , 2l^cos2^(m? — rrin) , ,\ , , 

^ ^ (9(71 + (72 + as + (74) + ^(12ai + 2(72 -2a3) h'~ 1.3 -10-2 



2P V X ■ - ■ o . . ^ 



(4 



Ignoring subleading MSW effects V, for simplicity, and considering oscillation lengths 
t = L ^^a^^a) ; we then observe that the independence of the length L result of the 
experimental fit, found above, may be interpreted, in this bounding the stochastic 

fiuctuations of the metric to Qui + (72 + c^s + ~ 1-3. ■ 10~^. This is too large to be a quan- 
tum gravity effect, which means that the contributions to the damping due to stochastic 
fluctuations of the metric, as in the model of |3| above, cannot be the explanation of the 
flt. 
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B. Conventional Explanation: Energy Uncertainties 



The reader's attention is called at this point to the fact that such time-reversal invariance 



decoherence may also be due to ordinary uncertainties sl in energies and/or oscillation 
lengths, which are unrelated to quantum gravity effects. For instance, consider the ordinary 
oscillation formula for neutrinos, with a mixing matrix [/, 

-Po/3 = Pal3{L, E) = 

n n / A 2 



^^^-^Y. H Re {u:,UpaUo.bU;,) sin' ' 

a=l I3=l,a<b 
n n 



AE 



2J2 E l^{KaUpaU.bU;,)sinl^^] (4.10) 

a=l b=l,a<b ^ ^ 

where a, P = e, fi,T, a,b = 1, 2, ...n, Am^j, = — m^. 

In general there are uncertainties in the energy E in the production of a z/ (and/or z7)- 
wave, and also in the oscillation length. As a result of these uncertainties one has to average 
the oscillation probability ()4.10|1 over the L/E dependence. 

Considering a Gaussian average jQl, 

/°° 1 i^-ef 

dx P{x)^^e-^^ 
■oo 0"v2vr 

^ = (x), a = a/ {{x — {x)Y, X = L/AE, and approximating {L/E) ~ (L) / (E) we obtain 

{Pais) = Saf3 — 



n n 

2E E Re (f/l^^,^„,f/;,) (l - cos(2Mm^,)e-2'^^(^--)^ 

a=l /3=l,a<fe 
n n 

-2E E Im (f/lf//5a^„.f/;,) sin(2Mm^,)e-2'^^(^--)^ (4.11) 

a=l b=l,a<b 

Notice the exponential damping factors due to the fluctuations a. 

In fact, as discussed in jsf, there are two kinds of bounds for a: A Pessimistic: 

one, according to which a ~ Ax ~ A^ < + and an Optimistic: a < 

,rx / \ 1/2 

4(i?> l^i(L>J + V(,E)\ J ■ 

In our case, where we consider long baseline experiments, the uncertainties in the oscil- 
lation length L are negligible, and hence the two cases degenerate to a single expression for 

^ _ {L) AE 
" ~ 4(E) E ■ 
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The damping exponent, then, in ()4.1ip . arising from the uncertainties in the energy of 
the (anti)neutrino beam, becomes 

2a\/\m^f = 2jj^ ( ^) ' {Amy . (4.12) 



(my V E 

As mentioned above, for oscillation lengths we have LArn^ /2E ~ 0{1), and hence, the result 
of the best fit ()3.3p . (13.41) . implying independence of the damping exponent on L (irrespective 
of the power of L), yields an uncertainty in energy of order 

AE , , ^ 

—^-1.6-10-^ (4.13) 

E 

if one assumes that this is the principal reason for the result of the fit. This is not an 
unreasonable result, thus implying that the result of the fit may be interpreted as being due 
to ordinary physics associated with uncertainties in the energy of the neutrino beam. 

The important difference from the stochastic fluctuations of gravity medium, discussed 
above, lies on the fact that the period of oscillation is not affected 1)4.111) by the above 
averaging procedure, in contrast to the stochastic gravity cases fl4.4|) and ()4.7|) . and thus in 
principle the effects can be disentangled. However, in general these latter corrections are 
small, and beyond the sensitivity of the current experiments. Nevertheless, as we have seen, 
some effects, such as the time-reversal symmetric stochastic fluctuations of the background 
metric ()4.7|) . can be already disentangled straightforwardly by their order as compared with 
the above energy-uncertainty effects due to ordinary physics. 

The precise energy and length dependence of the damping factors is an essential step 
in order to determine the microscopic origin of the induced decoherence and disentangle 
genuine new physics effects I2I from conventional effects, which as we have seen above may 
also contribute to decoherence-like damping jj]. For instance, as we discussed above, some 
genuine quantum-gravity effects, such as the stochastic fluctuations of the space time, are 
expected to increase in general with the energy of the probe j^, as a result of back re- 



action effects on space-time geometry, in contrast to ordinary-matter-induced 'fake' C 
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violation and 'decoherence-looking' effects, which decrease with the energy of the probe 
At present, the sensitivity of the experiments is not sufficient to unambiguously determine 
the microscopic origin of the decoherence effects, as we have seen above, but we think that 
in the near future, when experiments involving both higher energy and precision become 
available, one would be able to arrive at definite conclusions on this important issue. Thus, 
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phenomenological analyses like ours are of value and should be actively pursued, in our 
opinion, in the future, not only in neutrino physics but also in other sensitive probes of 
quantum mechanics, such as neutral mesons. 

V. CONCLUSIONS AND OUTLOOK 

In this work we have presented a complete analysis of three-generation neutrino transition 
probabilities, which include decoherence effects with guaranteed positivity. In our opinion 
this is the first complete, mathematically consistent, example of a Lindblad-decoherence 
model for neutrinos with full three generation mixing. 

We have shown that decoherence effects can account for all available neutrino data, 
including LSND results even in a minimalistic scenario (with only one new parameter, which 
parametrises all the decoherence effects). Contrary to other approaches in the literature, 

n nn 

using sterile neutrinos and following the spirit of our earlier work p, 112I|, we have 
attempted to interpret the LSND results not by means of oscillations, but as a decoherence 
effect inducing damping in the oscillatory terms, which is also present (as a per-mil-ish 
additional suppression) in other neutrino experiments as well. 

The specific oscillation length L dependence of the single decoherence parameter, im- 
plying L-independence of the corresponding damping exponents, could in principle find a 
natural explanation in some theoretical models of stochastic quantum gravity. However, its 
order of magnitude seems incompatible with this possibility, at least in the concrete space- 
time foam microscopic models considered here, since it would imply quantum-gravity effects 
unrealistically large. 

On the other hand, the result of our fit can find a natural explanation in terms of ordinary 
physics. It could be due, for instance, to uncertainties in the energy beam of (anti) neutrinos, 
and indeed this scenario seems to provide the most natural explanation of our fit. 

We now remark that quantum-gravity contributions could indeed be present, and lead to 
similar damping in oscillatory terms, but their suppressed order of magnitude would imply 
that they could only be probed at higher energies. For instance, high-energy neutrinos 
detected from distant supernovae, may probe these issues further, since they will increase 
significantly the sensitivity to genuine quantum gravity effects [2^, and thus may probe the 
induced changes in the damping as well as the oscillation period, as discussed in this work. 
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It goes without saying that there is much more work to be done, both theoretical and 
experimental, before definite conclusions are reached on this important issue, but we believe 
that neutrino (astro)physics will provide a very sensitive probe of new physics, including 
quantum gravity, in the not-so-distant future. 
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